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ON  BUSEMANN'S  DIAGRAMS  FOR  STEADY  (STATIONARY)  SHOCKS. 

1.  Linear  shock. 

Throughout  we  shall  assume  that  the  fluid  is  in  a 
steady  state.   Presently  we  assume  it  to  fill  a  straight 
tube  so  that  its  state  can  be  described  by  one  velocity 
component  u,  and  by  thr  pressure   p   and  the  density  p  . 
We  assume  that  pressure  p   and  density  P     are  connected 
through  a  relation 

(1.1)  p  =  <r(s)  f(P) 

s  denoting  the  entropy  of  the  state  (-;;-) .   The  sound  speed 
is  then  given  by 

(1.2)  c  =  [y  p/p  ]^/^  . 

We  also  introduce  the  enthalpy   i,  and  consider  it  also  a 


r  2 
-  \  c 


function  of  0   and   s   (•■>-), 

(1.3)        i  =  (r(s)  g((>)  =  j  -S-  ^p 

The  quantities  p  and   i   are  related  through 

(1.4)     P  ^l/<)p   =  ^p/(^P  =  c^,  i.e.   pg'(f^)  =  f'(P) 

If,  in  particular,  the  gas  is  ideal, 

(1.1)  P  =  cr  (s)y  ,   1  <  6''  , 
we  have 

(1.2)  c  =  [  iTP'^p]^/^ 
and 

(1.3)  1  =  -T-^-r  r^p  "  ^^/^  ^  -  1)  . 

Further  the  "Bernoulli  quantity" 


(-;;-)  We  assume   f(P)  >  0,  f'(^)  >  O",  f"(\^)  >  0,  flO)  =  0, 
f(oo)  =00;  (r(s)  >  0,  (r'(s)  >  0,  (r(-.'>=')  =  0,  cr(o<£')  =00  . 

(-i:-)  We  assume   g(0)  =  0;  then  <?g(p)  >  f(f )  follows  from  (1.4) 
Also  relation  (1.4)  implies  d   1/5'p  >  0. 
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(1.5)  ^  ^  I  ^^ 

is  introduced  and  all  states  are  then  considered  which  have 
the  sQme  entropy  s  =  s  and  the  same  "Bernoulli  constant" 
i°  : 

(1.6)  1  .1^2  =1^2  ^  1°   . 

When   3   and   i   arc  prescribed,  to  every  value   u  of 
the  velocity,  u  <  c  ,  the  quantities   P  and  p  are 
determined,  since  the  relation  (1.3)  can  be  solved  for  o 
by  virtue  of  (1.4).   In  particular   c  and  p   are  then 
definite  functions  of  u  : 

(1.7)  c  =  c(u)  =  c(u,  s)    (-"-) 

(1.8)  p  =  p(u)  =  p(u,s)    ih      . 
Especially,  for  ideal  gases, 

(1.7)  c  =  [ g (c   -  u  )]  ^ 

(1.8)  p  ^    (Ti-l/i  ^   -  1)]^  J^  -  1  (^2  _  ^2^^  d'/i^     -   1) 

Among  these  states  there  will  be  one  for  which 
the  sound  speed  coincides  with  its  velocity,  i.e.  for  which 
c(u)  =  u.   This  velocity  is  called  "critical  sound  speed". 
Evidently  we  have 

(1.9)  21°  =  Z^   =  -f-^-c"^   . 

If  u  <  c"  Vie   also  have  u  <  c(u)  and  the  flov;  is 
called  subsonic,  if  u  >  c"  we  also  have   u  >  c(u)   and  the 
flow  is  called  supersonic.   It  is  clear  that  it  will  bo 


(--■)  By  virtue  of  the  assumptions  in  footnote  (---), p.  1,  we 
h^ve  ^a^c/^^^^  >   0.  ,.-. 

(-;;-)  By  virtue  of  the  assumptions  in  footnote  (-"■),  p.l,  the 
function  p(u,  s)   decreases  v/he-n  u   increases. 
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convenient  to  use   c"   as  a  referonce  velocity. 

We  note  a  few  more  relations  needed  for  the 
discussion  of  the  function  p  =  p(u)  as  given  by  (l.S), 


(1.10) 


-  u 


Prom  (1.5)  we  have 

from  (1.4)  we  then  obtain 

(1.11)  Jr=:.j.uA2 

from  (1.2)  therefore 

and 

(1.13)        ^-^  - 


P  u 


-  -p  =  -  P(l  -  n/c). 
The  last  relation  shows  that   the  curves   p  =  p(u)  have 
negative  curvature  for  subsonic,  positive  curvature  for 
supersonic  velocity;  the  point  of  inflection  is  at  the 
critical  velocity.   It  is  to  be  observed  that   separation 
of  the  two  regions  u  >   c"   is  independent  of  the  entropy 


and  dependent  only  on   i  ., 


/ 


/ 


/ 


\ 


— -?■, 


-C 


c 


£ 


By  virtue  of  (1.12)  v^e  can  say  that  the  negative 
slope  of  the  curve   p  =  p(u)  represents  the  mass  crossing  ^ 
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unit  cross-section  per  unit  time.   A  significant  quantity  is 
also  the  intercept  on  the  p-axis  of  the  tangent  line  at  a 
point  of  the  curve  p  =  p(u).   We  find  for  this  intercept 

p  +  ^u   , 
i.e.  the  sum  of  the  pressure  and  momentum  crossing  a  unit 


area  per  unit  time. 


It  should  further  ho  mentioned  that  the  functions 


U. 


p  =  p(u)  for  different  entropies  differ  only  by  a  constant 
factor. 

It  is  now  very  easy  to  represent  a  shock  transi- 
tion graphically.   We  distinguish  the  left  or  the  right     ' 
hand  side  of  the  stationary  shock  by  (/)  and  (r).   Assume 
that  (/)  be  the  front  side  of  the  shock  and  that  uy  >  0. 
The  shock  relations  are 

(1.14)  p^u^  =   ^^u^ 

(1.15)  p^  +  ^^u|  =  p^  +  J-'^u^    , 

while  the  third  relation  can  simply  be  written  as 

(1.16)  i^  =  iy   or  c^  =  Cy 

The  two  points  in  the   (p,u) -plane  corresponding  to 
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both  sides  of  the  shock  may  be  located  on  tv/o  curves  p  =  p(u) 

having  the  same  value   c"   but  differing  in  entropy. 

The  relation  (1.14)  can  now  be  interpreted  to  say 

that  the  tangent  lines  at  the  two  p  =  p(u)  curves  at  the 

points  (/)  end  (r)  have  the  same  slope,  and  relation  (1.15) 

shows,  in  view  of  what  was  said  just  before,  that  those  two 

tangent  lines  have  the  same  intercept  on  the   p-axis.   That 

is,  the  two  tangent  lines  coincide;  or,  the  two  points  (X) 

and  (r)  lie  on  a  common  tangent  line  for  the  two  curves 

p  =  p(u)  v;ith  the  same  value  of   c".   It  is  now  immediately 

seen  by  use  of  the  curvature  properties  of  the  curve 

p  =  p(u)  that  to  every  point  iX)    with  uy  >  c"   there  is 

just  one  point  (r)  having  the  same  tangent  line,  clearly 

u  <  c".   Also  obviously 
r  ** 


which  implies 
(1.17)    f 


p^(u)  >  py(u)    , 


S    <  Sy 


\ 


w 


We  mention  that  for  ideal  gases  the  tv;o  velocities 


u  and  Uy  are  related  through 
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-;:-2 

(1.18)  ^/•'^r  ~  ^  ' 

this  follows  immediately  when  one  divides  relation  (1.15)  by 
(1.14)  and  eliminates   P'-'-p   through  (1.3),  (1.6),  (1.9), 
and  (1.16)  . 

Another  graphical  representation,  which  perhaps 
may  be  useful,  is  obtained  by  Legendre  transformation.   One 
introduces 

(1.19)  q  =  P  +  Pu^ 

as  a  function  of  the  Independent  variable  pu.  As  is  seen 
from  (1.13),  this  is  possible  for  u  <  c"  and  for  u  >  c". 
The  domain  of  Pu  is  then 

(1.20)  0  <   Pu  <  0'"'c""'   , 
where  P"  =  P(c")  . 

We  obtain  two  branches  for'   q  =  q(  P  u),  a  "lower" 
branch,  corresponding  to  u  >  c",   with  q(0)  =  0,  and  an 
"upper"  branch,  corres  ponding  to  u  <  c"  ,  with  q(0)  =  p  , 
For  different  entropies  the  functions  q  =  q(Pu)  differ  only 
by  a  constant  factor.   The  lower  branches  for  different 
entropies  do  not  intersect  each  other,  and  the  same  holds 
for  the  upper  branches  for  different  entropies.   However,  it 
is  immediately  seen  that  through  every  point  on  a  lower 
branch  one  can  pass  an  upper  branch  belonging  to  a  higher 
entorpy.   Such  an  intersection  of  lower  and  upper  branches 
for  different  entropie^  corresponds  to  n  shock  transition, 
as  is  o  vious  from  relations  (1.14)  and  (1.15). 
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2.  Two-dimensional  shock 


))■,,  » 


-/»*1>L 


pU- 


f*lL*       /°<^ 


1°  =  I  c^  =  i  +  l(u^  +  v^) 


To  Investigate  two-dimensional  "oblique"  shocks 
one  may  also  consider  the  manifold  of  all  states,  given  by 
velocities  (u,v),  by  p   and   P  or   s,  which  have  the  same 
entropy   s  and  the  same  Bernoulli  constant 
(2.1) 

c  and  p  are  then  functions  of  u   and  v,  as  a  matter  of 

/   2    2 
fact,  functions  of  /  u  +  v   ,  which  reduce  for  v  =  0  to 

the  functions  of  u   as  considered  in  the  preceding  section. 

In  the  (u,v) -plane  we  consider  the  circles 

2  ^   2  ^  •;;-2  consisting  of  subsonic  states  and  the  ring 
u  +  v  <  c  ^ 

,;;.0       2      2       2 

c"   <  u  +  V  <  c   containing  supersonic  states.   Also  we 


observe  as  before 


(2.2) 
(2.3) 


£  =  _ 
u 


? 


u. 


v 


-  fv, 


^2„      „,.     2,2.  ^  2^  .  .  ^  „  .  ,2/,e,^ 


_P   _ 


b    U 


-  j3(i  .  uVc  )   , 

O     p  _  uv 
gUTv    ^2 


a 


V 
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y'o>.   /<)vv  -  .2,.   u2 


2 


Thus  we  see  that  the  curvature  is  positive  or  negative  for 
subsonic  or  supersonic  velocities,  respectively.   The  sub- 
sonic and  supersonic  parts  of  the  surface   p  =  p(u,v)  may  be 
called  the  "dome"  ai  d  the  "rim"  of  the  surface. 

Suppose  now  two  states  {/)    and  (r)  are  separated 
by  a  stationary  straight  line,  an  "oblique"  shock.   Let  u 
and  u.  be  the  normal  and  the  tangential  components  of  the 
velocitv.   Then  the  transition  relations  are 


n 


(2, 

5) 

^Xt 

=  ^rt 

(2. 

6) 

fxVn 

=     pu 
r   rn 

(2. 

7) 

P/ 

^fX^Xn 

^r 

and  further 

(2. 

8) 

^5 

=   ^r   = 

.  0 

1        or 

fr- 


rn 


'^r 


Relation  (2.':.)  says  that  the  vector  (uy-u^,vy-v^) 
is  normal  to  the  shock  direction,  i.e,  the  shock  direction 
can  be  found  as  the  normal  to  the  line  connecting  (uy, vy) 
and  (u  ,v  )  in  the  (u,v) -plane.   In  the  (u,v,p) -space  the 
two  points  (X)  and  (r)  are  located  on  two  surfaces  p  =  p(u,v) 
having  the  same   c'",  differing  only  in  entropy. 

It  is  now  easily  seen  that  the  straight  line  in  the 

(u,v,p)-space  connecting  the  two  points  uy,vy,py  and 

u  ,v  ,p   is  a  common  tangent  to  the  two  surfaces   p  =  p(u,v). 

Indeed,  the  relation 

\ 
(2.9)  PX  -  Pr  -^  X^ 


""  (-X  --r^'y^ 


X 


(Vy  -  V^)   =  0 
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is  an  immediate  consequence  of  -y-^   =  -  p  u,  -r~^   =  -  p  v,  and. 
the  three  shock  transition  relations. 

We  consider  now  one  state  (0)  and  ask  for  all  states 
(l)  that  may  be  connected  with  (0)  by  an  oblique  shock.   The 
straight  lines  in  the  (u, v,p)-space  connecting  (0)  v/ith  all 
such  states  (1)  will  then  lie  on  the  tangent  plane   T   to 
p  =  p  (u,v)   at  the  point  (O).   The  procedure  is  now  clear, 
Vi/e  consider  a  surface   p  =p(u,v)with  any  entropy  and  inter- 
sect it  with  the  tangent  plane   T  .   In  general  there  v;ill  be 
two  tangent  lines  at  the  intersection  curve  passing  through 
(O),   The  points  of  contact  of  these  tangents  with  the 
intersection  curve  are  states  (l)  that  can  be  connected  with 
(0)  by  a  shock. 

To  investigate  in  detail  the  locus  of  the  possible 

points  (1)  we  first  assume   u  >  c",  i.e,  we  first  assume  the 

state  u   to  be  supersonic.  We  then  consider  the  surfaces 
o 

p  =p(ujv)for  decreasing  entropy.   When  the  entropy  is 

infinite,  p  =  0;  when  the  entropy  decreases  the  dome  of  the 

surface   p  =  p(u,v)  will  rise,  touxsh,  and  eventually  cross 

the  tangent  plane   T   (-:>),   The  state  at  which  the  dome 

touchou"  (1?    corresponds  to  the  straight  shock  as  considered 

in  the  preceding  section.   When  the  dome  has  crossed  the 

plane   T   the  intersection  will  be  an  oval  and  clearly  there 
o 

will  be  two  tangent  lines  yielding  two  states  (1)  that 

correspond  to  tv/o  possible  oblique  shocks.   When  the 

entropy  decreases  to  s  ,  these  two  points  move  into  the  point 

(0),  So  far  the  state  (0)  was  on  the  front  side,  the  state 

(-;:-)  It  can  be  seen  that  at  the  intersection  of  the  rim  with  the 
plane  T   there  are  no  tangent  lines  through  (O)  as  long  as 
the  entropy  is  greater  than  s   * 


lOi 


(l)    Qt   the      back   side   of   the    shock.      V.Tien      s     has  become    less 

than     a.,   this   situation   is   reversed.      The    state    (1)    v/ill 

2    2    2 

then  have  greater  speed,   u,  +  v,  >  u   ,  and  eventually 

approach  the  Intersection  of  T  with  the  limit  circle, 

2    2     2 
i.Ci  u^  +  v^  ->  c  .   A  similar  discussion  could  be  carried 

out  for  u   <  c" :  the  state  (0)  is  then  of  course  alwavs 
o 

on  the  back  side  of  the  shock. 

The  locus  of  all  points  (1)  which  may  be  connected 
with  (0)  by  a  shock  -  or  their  projections  in  the  (u,v)-plane" 
is  called  shock-polar.   In  case  of  an  ideal  gas  its  equation 
is  easily  found.   Let  the  shock  direction  be  given  by  the 
angle  /3  with  respect  to  the   u-axls.   Then  we  have 
(2.10)  Uq^  =  Uq  cos^  , 


(2.11) 


u   =  u   sln|3 
on     o     r 


^ 


—^ 


p  C{ron-t-i 


^0 


I  I  I  i  '  '  I   '  /  /  I  !  '  I   '  '  '  I  '  1 1  ■ 


^ 


-^■if- 


h^ 


\^^ 


1  \ 


-VA-J 


r<^ 


For  the  state  (1),  omitting  the  subscript  (l),  v;e  have  from 

(2.5),  (2.6) 

(2.12)  u^  =  Uq^  =  u^  cos^  . 

Dividing  (2.7)  by  (2.6)   and  eliminating  p"  p  by  (1.3)   , 

(2.1)  ana  (2,8)  one  obtains 
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n  on 

°^  -2  _  ^2  .2    2 

(  'J'  -  1)~~—  +  (  ^  +  l)u^=  (a^  -l)L^-^Qt+  (3^  +l)u^^ 
n  on 

We  multiply  by  ^^-^^^  ,  divide  by  u^  -  u^^  , 

2    2 
ma|clng  use  of  u.  =  u  .  .      Then  we  obtain 

(2.13)  u  u  =  c"^  -  -T-r-T  ^r.^-     > 

on  n         j^  +  1  ot   ' 

a  relation  which  is  m alogous  to  (1.18).   From 
u  =  u,  cos/3   +  u  sin|3 

V  =  -u.sin^   +  u   cosP 

p.nd    (2.10),     (2.11),     (2.12)    we    obtain   tho   pr,rr.metric 
representation 

(2.14)  u  =  -^  u^    cos^/j     ^  §- 

V  +   [{  /^  ^   cos^p      -   1)Uq  +   c''''Vuo]    cot^    . 

Elimination  of  f>   yields  first 

(2.15)  V  =  (u  -  u^)  cot  p 

an  d  then 

/OT/'\  2_/  n2u-c       /Uo 

(2.16)  V  =  (u  -  u  )  — ijj riT^ 

oo 

with 

2u      -;;-2 

The  locus  of  this  equation  is  a  strophoid  as  indicated  in 
the  diagram. 

Busemann  considers  anly  the  "backward"  bcanch 
for  u  >  c\    i.e.  the  branch  consisting  of  those  states 
(1)  which  may  be  connected  with  (O)  in  such  a  way  that  (O) 
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—  ^ 
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CP        !'■ 


is  on  the  front  side,  (1)   on  the  back  side  of  the  shock. 
The  other  branches  should  however  not  bo  forgotten. 

We  observe:  Fnen  u   approaches  ■,  p",  the  backxvnrd 
branch  shrinks  to  the  point  (c'"',0).   VVhen  u^  approaches 
c,  the  value  u   also  approaches  c;  the  backivard  branch 
approaches  the  circle   v  =  (c  -  u)(u  -  c"  /c ) ,  while  the 
forward  branches  approach  the  line  u  =  c. 


13. 

The  bactavp.rd  branch  of  the  shock  polar  may  easily 
be  used  to  determine  the  oblique  shock  v;hen  the  direction  of 
velocity  on  the  back  side  is  given  by  the  angle  od  with  the 
u-axis.   One  need  only  draw  a  line  in  that  direction  through 
the  origin  in  the  (u,v) -plane  and  intersect  it  with  the 
backward  shock  polar.   There  will  be  either  two  intersec- 
tions, one  intersection,  or  none  at  all.   If  the  point  of 
intersection  is  found  on  the  polar,  the  shock  direction  is 
found  as  the  perpendicular  to  the  line  connecting  (0)  and 
(1).   The  normal  and  tangential  components  of  both  velocities 
are  also  easily  read  off  from  the  diagram. 

It  is  interesting  t6  investigate  the  limit 

position  oC  of  the  angle  oL    .   This  angle  depends  on  the 

initial  state  (O);  if  the  initial  velocity  u   approaches 

c",  the  angle  oC       approaches  zero.   If,  however,  the 

velocity  u   approaches  c,  the  angle  o^     approaches  an 

"ultimate"  value  given  by 

sin  oL^  =  1/  ^    ; 

in  particular 

o^    =  45.5°      if     '3^  =  1.4      . 
o 

The  corresponding  limit  angle   /^   for  the  shock  direction 
is  given  by 

^  o   4    2 
The  "ultimate"  limit  value  j5  is  therefore 

/3  =  67.75^ 
'  o 

It  should  also  be  noted  that  the  state  correspon- 
ding to  the  limit  angle  is  subsonic. 

The  angle  oC  under  which  the  polar  intersects  the 
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2  2  2  ^ 

circle     u     +  v     =   c       may  be  denoted  by    oL   ;    this   r.ngle 


approaches   0  when     u 


when  u  — >  c 
o 


c  or  0  and  an  ultimate  value  ^ 


Further:  If  cx^  ^»  0,  the  state  u  approaches 
either  to  that  corresponding  to  the  straight  shock  or  to  the 
state  (0)  itself.   In  the  latter  case  the  shock  direction 
approaches  the  characteristic  direction  vrhich  is  therefore 
perpendicular  to  the  tangents  at  the  double  point. 

It  is  of  course  possible  to  apply,  as  in  one- 
dimensional  case,  a  Legendre  transformation  and  to  consloer 
the  two  bran cos  of 

q  =  p  +  P (u^  +  v^) 
as  functions  of  p  u  and  Pv;  one  then  has 

^  q/<)  (p  u)  =  u,  ^  q/d  (  Pv)  =  V. 


Again  the  line  connecting  two  points  that  correspond  to  both 
sides  of  a  shock  is  tangent  to  both  surfaces   q  =  q(Pu,Pv), 
The  shock  polar  connecting  a  supersonic  state  with  (1)  super- 
sonic (ii)  subsonic  states  is  indicated  in  the  diagram. 
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3,  Remarks  on  shock  reflec^Torfr 


Let  us  consider  the  steady  state  flow  along  a  v/all 
crossing  two  shocks,  as  indicated  in  the  diagram. 


To  investigate  the  possibility  of  such  a  flow  one  may 

consider  (u  ,v  )  and  P  >P    £^s  given,  rssuming  that  ('^q»'^q) 

is  supersonic.   Then  one  need  only  draw  the  shock  polar  for 

the  point  (u  ,v  )  and  observe  whether  this  polar  intersects 
^       o'  o  ^ 

the   u-axis  in  its  forward  and  its  baclovard  branch.   It  is 


16. 


clear  that  an  Intersection  v/ith  the  backward  branch  will  exist 

only  if  V  /u  <  tan  o^  ,  v^fhere  o^^    is  the  limit  angle 
•'       0   0.       o  o  ° 

belonging  to  the  state  (O),   The  intersection  with  the 
forward  branch  will  exist  only  if  od  <  al      (see  p. 13), 
4.  Remark  on  the  Mach  effect. 

Consider  the  steady  flow  as  indicated  in  the 
diagram;  the  states  (1)  and  (2-)  are  connected  by  an  oblique 
shock.   On  the  other  hand  state  (1)  passes  through  an 
oblique  shock  over  to  a  state  (0)  and  (O)  passes  over  to 
a  state  (2)  through  an  oblique  shock. 


(t) 


m 


We  shall  show  that  such  a  cd>nf iguration  is 
impossible  ,   Consider  the  points  (O),  (1),  (2)  in  the 
(u, v,p)-space  and  the  surfaces  p  =  p  (u,v),  p  =  p^(u,v), 
p  =  p  (u,v)  belonging  to  them.   The  lines  connecting  points 
(1)  and  (2),  (1)  and  (0),  (0)  and  (2)  are,  as  was  shown, 
tangent  to  those  surfaces,  respectively.   The  plane  spanned 
by  these  three  points  is  therefore  a  tangent  plane  to  all 
three  surfaces.  Clearly  the  intersection  of  a  tangent 
plane  with  the  (u,v) -plane  is  per5e*ndicular  to  the  vector 
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